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ABSTRACT 

We compare the environmental and star formation properties of far-infrared detected and non- 
far-infrared detected galaxies out to z ~ 0.5. Using optical spectroscopy and photometry 
from the Galaxy And Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS), with 
far-infrared observations from the Herschel-AILAS Science Demonstration Phase (SDP), we 
apply the technique of Voronoi Tessellations to analyse the environmental densities of individ- 
ual galaxies. Applying statistical analyses to colour, r— band magnitude and redshift-matched 
samples, we show there to be a significant difference at the 3.5a level between the normalized 
environmental densities of these two populations. This is such that infrared emission (a tracer 
of star formation activity) favours underdense regions compared to those inhabited by exclu- 
sively optically observed galaxies selected to be of the same r— band magnitude, colour and 
redshift. Thus more highly star-forming galaxies are found to reside in the most underdense 
environments, confirming previous studies that have proposed such a correlation. However, 
the degeneracy between redshift and far-infrared luminosity in our flux-density limited sam- 
ple means that we are unable to make a stronger statement in this respect. We then apply 
our method to synthetic light cones generated from semi-analytic models, finding that over 
the whole redshift distribution the same correlations between star-formation rate and environ- 
mental density are found. 

Key words: method: data analysis - galaxies: statistics - galaxies: submillimetre- galaxies: 
star formation. 
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1 INTRODUCTION 

1.1 The environmental influence on galaxy formation and 
evolution 

It is clear that if we are to understand the process by which galaxies 
form and evolve, we have to consider the role that their immedi- 
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ate environment plays. 'Dres slerj ( p^80| was the first to show that 
there is a correlation between the morphology of a galaxy pop- 
ulation and the density of its environment. Further studies have 
since shown that disk dominated 'late-type' galaxy morphologies 
with high star formation rates (SFR) dominate underdense regions, 
while their elliptical 'early-type' counterparts, with low SFR, dom- 
inate the densest regions (^Postman & Geller^l98 4 , Dressie r et al. 
™^ jPommguez et al. 12001 [[Goto et al |2003[|K auffma nn et al.| 



emission ( Schmitt et al. 2006 , da Cunha et al. 2008 ; Nardini et al. 
2008, Bendo et al. 2010 2012, Groves et al. 2012, Smith eFaL 



1997 



2004 



O'Mill et al.|2008l|Lee et al.|2010[[^^jesinghe et al.|20T2l ). 



It has also been shown, with the advent of large area sur- 
veys such as the Sloan Digital Sky Survey (SDSS; |York et aT] 
|2000j ), that these galaxies can be categorised into two distinct op- 
tical colour populations, 'Red' and 'Blue\ where the colour of a 
galaxy is dependent on several internal properties that represent the 
evolutionary history of the galaxy; metallicity (Z), star formation 
history (SFH) and dust attenuation {A). These two colour popula- 
tions show, at a fixed luminosity, a correlation with density such 
that the densest regions are populated by the red, early-type pas- 
sive galaxies, with blue, star-forming late-types observed in less 
dense regions ([Poggianti et al.| [2006 ). However, earlier work by 
[Balogh et al.| ( |1997||1998| ) compared galaxies with similar lumi- 
nosities and morphologies from both dense cluster and low den- 
sity field environments, and found that that the SFR was still lower 
in dense cluster regions and thus the SFR-density correlation still 
held regardless of the morphology of the galaxy. This indicates that 
the observed SFR-density relation cannot be exclusively tied to the 
morphology-density relation; other processes must be influencing 
the observed correlations. It is currently believed that this reduc- 
tion in SFR with increased environmental density is directly linked 
to the stripping of cold gas from galaxies via some type of direct 
interaction, and several mechanisms have been invoked to explain 
this observed correlation. 

For example, major mergers ( [Barnes & Hernquist|1992| ) can 
cause a burst of star-formation activity and feedback from such star- 
burst events is able to prevent gas cooling, and as a result the gas re- 
mains out of pressure equilibrium with its environment. Due to this 
pressure difference the gas expands out of the central regions of the 
galaxy sweeping up the inter- stellar medium (ISM). This ejection 
of the ISM from the merger remnant can lead to further suppression 
of star formation ( jMac Low & Ferrara|19"99l|Gay et ar|2010j . Al- 
ternatively, other processes such as harassment, strangulation and 
ram-pressure stripping may also play an important role (see e.g. 
[Boselli & Gavazzi,2006, for a review). 

1.2 Far-infrared emission as a tracer of star formation 

Star formation within a galaxy typically increases the dust content 
of the ISM through processes associated with the short-lived mas- 
sive stars that inhabit these regions, such as supemovae, that re- 



distribute material into the surrounding ISM ( Dunne et al.||2003[ 
[Sugerman et al.|2006[ [Dunne et al.|2009l [Gomez et al.|2012| ). This 
dust then absorbs a significant fraction of the ultra-violet (UV) light 
emitted by the young O-B type stars associated with these regions 
and is heated to temperatures of around 20-40 K, emitting thermal 
radiation at far-infrared wavelengths. This makes the use of far- 
infrared emission from a galaxy a widely used diagnostic for the 
obscured SFR of a galaxy (K ennicutt 1998, Hirashita et a l. 12003 [ 
Driver et al.|[2007l[Cort ese et al. 2008, Nordon et al.||2010l]Buat| 



etal.|2010|[Dunne et al . 201 1 ; Smith et al. 2012al. 



However, other contributions to the UV radiation field which 
heats the dust, such as AGN and older stellar populations within 
the galaxy, may lead to overestimates of the SFR using far-infrared 



[2012a|[Smith et al.|2012"bl ). Conversely, in galaxies where the ISM 
is optically thin at UV wavelengths, the measured SFR will emerge 
directly from the UV and not in the far-infrared. In these galaxies, 
deriving the star formation from far-infrared emission may lead to 
an underestimation of the total SFR ( Kennicutt 199^). However, 
as more than 50 per cent of energy ever radiated from stars has 
been absorbed by dust and re-radiated into the infrared ([Puget et aL 



1996) [Fixsen et al.1p^98l [Adelberger & Steidel[[2000| [Finn et a" 



2010 1, with the bulk of star formation since z — 1 occurring in 



dust obscured galaxies ( 


Calzetti & Heckman|1999|[Le Floc'h et aL[ 


[20051 [Patel et al.[[2013 


1, only AGN, low metallicity systems and 



very passive but dusty galaxies will lie off the far-infrared to SFR 
relation. 

Initial studies of the relationship between SFR and infrared 
emission from galaxies focused on shorter wavelengths using the 
Infrared Astronomical Satellite (IRAS; [Neugebauer et al.|1984| and 
more recently the Spitzer Space Observatory ( fRieke et al.[[2004| . 
IRAS surveyed the vast majority of the sky between 12 — 100 /xm, 
providing a large census of dusty galaxies in the local Universe. 
Using these data, [Goto| ( [200 5 ) investigated the optical properties of 
4248 infrared- selected galaxies by positionally matching data from 
the IRAS with optical data from the SDSS. Using a volume limited 
sample at 2; < 0.06 and applying a 5th-nearest neighbour den- 
sity estimate, a trend was found such that galaxies with the highest 
infrared luminosities reside in relatively low-density local environ- 
ments, suggesting that star-forming galaxies favour underdense re- 
gions, in agreement with previous studies at other wavelengths. 

The environmental densities of //?A5-detected luminous infra- 
red galaxies (LIRGs; 10^^ < Lfir < 10^^ L©) at 0.03 < z < 
0.17 were also studied by Hwang et al.[( [2010| ). They found that the 
fraction of LIRGs was strongly dependent on both the morphology 
and the distance to the nearest neighbour galaxy. They conclude 
that the evolution of the SFR-density relation from high to low red- 
shifts is consistent with the idea that galaxy-galaxy interactions and 
merging play a critical role in triggering star formation in LIRGs. 

Additionally, Tekola et al. (2012 ) examined the relationship 
between star formation and the environments of LIRGs selected 
from IRAS and compared these with other types of high- and low- 
redshift galaxies out to 2; ^ 1. They identified that there was a 
luminosity (Ljr ^ 10^^h~^ L©) at which infrared selected galax- 
ies preferentially resided in higher density environments, compared 
to "normal" galaxies. Above this luminosity the average density 
increases, whereas below this luminosity, infrared- selected galax- 
ies reside in environments of equal density, similar to the general 
population. They conclude, therefore, that infrared activity for non- 
LIRGs is not dependent on density and that the SFR-density re- 
lationship for these galaxies is similar to that of blue galaxies at 
z - 1. 

At higher redshifts, |Feruglio et al.[ ( [2010l ) used 24/im observa- 
tions from Spitzer to investigate the environmental effects on star 
formation in LIRGs and ultra- luminous infrared galaxies ( ULIRGs; 
> 10^^ Lo) in the Cosmic Evolution Survey (COSMOS ;|Scoville 
[et al.|2007| ) at 0.3 < z < 1.2. They found the fraction of these 
galaxies to decrease with density out to z ^ 1, but that the relation- 
ship flattens out with increasing redshift. 

Due to the wavelength coverage of IRAS (12 — lOO/xm), the 
majority of galaxies detected by these studies were found by [Breg-[ 
[man et aL] (pr998 ) to be spirals and starbursts in the local universe 
(z < 0.1). This restriction resulted in the IRAS providing little in- 
formation about the cooler dust, which traces the bulk of the dust 



© 2012 RAS, MNRAS 000, ??-?? 



The environmental density of far-infrared bright galaxies 3 



mass (e.g. |Dunne et al.|2 01l1, in other galaxy populations, espe- 
cially early-type morphologies. In comparison, Spitzer can probe 
longer wavelengths (24 — IGO/xm) and therefore is less susceptible 
to this bias, although galaxies with the coldest dust temperatures 
would still be missed (Bales et aL]|2010[ [Symeonidis et al.||2011| 
^13 ). Considering that cold dust is present across all types of 
galaxies and is a major contributor to infrared luminosity ( Willmer 
|et al.|2009] ), and closely traces the total dust mass, it is crucial that 
we are able to select galaxies at longer wavelengths. 

With the launch of the Herschel Space Observatory ( |Pilbratt| 
|etal.|2010| ) we are now able to select galaxies at these longer wave- 
lengths. A number of studies have begun to investigate how star 
formation in a galaxy, traced by far-infrared emission at > 250//m, 
is linked to the environment in which the galaxy resides. |Dar-| 
|iush et all (|2011)) used far-infrared data from the Herschel As- 
trophysical Terahertz Large Area Survey (H-ATLAS; |Eales et al.| 
|2010| ) Science Demonstration Phase (SDP) to examine the ultravi- 
olet and optical properties and environments of low redshift galax- 
ies (0.02 < z < 0.2) from the SDSS and the Galaxy And Mass 
Assembly survey (GAMA; ^Driver et al.||20TT] [Hill et al.||20n] 
[Baldry et al.][20T0l ). They found that H-ATLAS detects predomi- 
nantly blue/star-forming galaxies, with a minor contribution from 
red galaxies (comprising highly obscured and passive systems). Us- 
ing the 5th-nearest neighbour as an estimate of the environmental 
density, they found that the fraction of H-ATLAS detected galaxies 
is much higher (^ 70 per cent) in low-density environments com- 
pared to high-density environments, where the fraction was found 
to be ^ 30 per cent. However, the detection rate of red and blue 
galaxies appears to be similar for both high- and low-density en- 
vironments, indicating that it is the colour of a galaxy, rather than 
the density of its local environment, that governs whether it is de- 
tectable by H-ATLAS. 

A consistent result was also found by [Rowlands et al.| ( |2012| ), 
who found that H-ATLAS detected early-type galaxies tend to have 
bluer (NUV-r) colours, higher SSFRs and younger stellar popula- 
tions than optically observed early -type morphologies. They com- 
pare 354 spiral and 30 early-type galaxy morphologies at low red- 
shift {z < 0.18), finding no significant difference between the en- 
vironmental densities of these populations. However, it is possible 
that they are not sampling a large enough range of environments 
with such small population samples. 

[Coppin et al.| ([2011) also used far-infrared data from H- 
ATLAS to examine the centres of 66 optically selected galaxy clus- 
ters at z ^ 0.25 to search for statistical evidence of obscured 
star formation in the cluster population. Using Voronoi Tessella- 
tions (described in Section |3.1| l to locate cluster members, they 
found an excess in the surface density of far-infrared sources within 
^ 1.2Mpc of the centre of these clusters. They conclude that the 
far-infrared emission is associated with dust-obscured star forma- 
tion in cluster galaxies, translating to a rate of ^ 7M0 yr~^. This 
SFR, maintained over the 3 Gyr since z — 0.25, would contribute 
enough mass to construct a typical SO- type bulge that would match 
the observed increase in bulge-dominated galaxies in cores of clus- 
ters over the same timescale. 

The effects of environment on the far-infrared properties of 
galaxies are also discussed by |Davies et al.| ( [2010) , who use data 
from the Herschel Virgo Cluster Survey (HeViCS; [Davies et al.| 
|2012f , finding that relatively few faint far-infrared sources that can 
be associated with confirmed Virgo cluster members. Furthermore, 
studies by'Cortese et al."("2010a"b) present Herschel observations 
of the perturbed galaxy NGC 4438 in the Virgo cluster and identify 
regions of extra-planar dust up to ^ 4 — 5 kpc away from the galaxy 



disk. This dust is found to closely follow the distribution of stripped 
atomic and molecular hydrogen, supporting the idea that gas and 
dust are perturbed in a similar way within the cluster environment. 

In contrast to these results, |Geach et al.| (|2011a'), using 
24/im observations from Spitzer, investigated large-scale filamen- 
tary structure surrounding rich clusters out to ^ ^ 0.55, and found 
that the SFRs of individual galaxy members within a cluster are 
not significantly different to identically selected field galaxies. Al- 
though pockets of enhanced star formation were observed, they 
suggest that this is the result of some 'pre-processing' effect where 
satellite groups have star formation triggered via gravitational tidal 
interactions during cluster infall. However, they state that there is 
no environmental mechanism acting to enhance the star formation 
within individual galaxies. 

It is evident that the majority of these studies have either 
used density measures that do not detect differences on the small- 
est scales (i.e. nth-nearest neighbour or aperture gridding) and/or 
they have focused entirely on narrow and local redshifts (z < 0.2). 
In this paper we use data from H-ATLAS to investigate the envi- 
ronmental dependence of far-infrared emission using a technique 
based on Voronoi Tessellations. Unlike the nth-nearest neighbour 
technique, Voronoi Tessellations calculate the environmental den- 
sity of galaxies on individual galaxy scales and hence can probe the 
environmental density to a greater degree of accuracy. 

In Section [2| we outline the optical and infrared data that we 
use. In Section |3] we present how both the spectroscopic and pho- 
tometric redshifts for our sample of galaxies were measured and 
sampled and introduce our algorithm to estimate the environmental 
density. In Section |4] we present the results of our analysis to de- 
termine whether there are any differences in environmental density 
between the far-infrared bright and faint sources, and investigate 
whether the SFR is linked to the environmental density. In Sec- 
tion |5] we compare our results to semi-analytic models and discuss 
the physical mechanisms that may explain our results. In Section|6] 
we discuss our results in the context of the physical mechanisms 
outlined above and in Section [7] we summarise our findings. We 
adopt a cosmology throughout with Qm = 0.30, Ha = 0.70 and 
Ho = 71kms"^ Mpc"\ 



2 OBSERVATIONS 
2.1 Far-infrared data 

We use far-infrared data from the science demonstration phase of 
H-ATLAS ( [Rigby et al.l[20TTl ). H-ATLAS provides data across a 
wavelength range of 100-500 /xm using the Photo-detector Array 
Camera and Spectrometer (R\CS; Pogli tsch et al. 2010) at 100 
and 160 /xm; and the Spectral and Photometric Imaging REceiver 
(SPIRE; |Gnffin et al.|2010 | at 250, 350 and 500/im. The H-ATLAS 
observations consist of two scans in parallel mode reaching 5a 
point source sensitivities of 132, 126, 32, 36 and 45 mJy in the 100, 
160, 250, 350, and 500/im channels, respectively, with beam sizes 
of approximately 9, 13, 18, 25 and 35 arcsec in the same five bands. 
The SPIRE and PACS map-making procedures are described by 
Pascale et al7(20ir) and Ibar et al.'('2010) respectively, while the 
catalogues are described by Rigby et al. (201 1 p . [Smith et al. ( 201 Ij ) 
used a likelihood ratio (LR) method to associate optical counter- 
parts with the H-ATLAS galaxies down to a limiting magnitude of 
r = 22.4 within a 10 arcsec radius. This resulted in optical counter- 
parts for 2, 423 objects from the H-ATLAS 250/im catalogue, each 
with a reliability R > 0.8 which ensures not only that the contam- 
ination rate is low but also that only one r— band source dominates 
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the far-infrared emission. While the entire H- ATLAS survey aims 
to compile a catalogue of ^-^10^ extra-galactic far-infrared sources 
out to z 3, the SDP field covered ^ 3 per cent of this, over an area 
of - 4.0 deg X 3.6 deg centred on (09^05"", -hO°30'). To maintain 
consistency with our Optical catalogue outlined in Section we 
limit our far-infrared catalogue magnitude to r = 21.5. 

2.2 Optical and near-infrared data 

We use spectroscopic redshifts from both the SDSS and the GAMA 
survey Data Release One (DRl). Spectroscopic redshifts are pro- 
vided for magnitude limits of r < 19.4 , K < 17.6 and z < 18.2 in 
the GAMA 9hr (G09) field which includes the H- ATLAS SDR This 
is combined with photometric redshifts derived from the combina- 
tion of optical (ugriz) SDSS and near-infrared (YJHK) UKIDSS- 
LAS imaging data as detailed in jSmith et al.| ( |20lT] ). This complete 
optical-near-infrared catalogue, containing photometric and spec- 
troscopic redshifts (hereafter named the Optical-9hr catalogue), 
totals 909,985 objects from which we remove all sources with 
r > 21.5 due to the fact that at fainter magnitudes the signal-to- 
noise ratio deceases to an extent where errors associated with the 
photometry become large. In addition, we remove objects classified 
as point-like in the SDSS imaging. This reduced the Optical-9hr 
catalogue to 323, 969 objects, of which 8, 875 had spectroscopic 
redshifts across a redshift range of < 2; < 1.2. 



3 ENVIRONMENTAL DENSITY MEASUREMENT 

In this section we describe our method of determining the environ- 
mental density of individual galaxies in redshift slices. First, as the 
vast majority of the galaxies in our sample do not have spectro- 
scopic redshifts, we are forced to use their photometric redshifts 
in order to establish where they reside in three dimensional space. 
Spectroscopic redshifts have errors of the order of Az ^ 10~^ 
( [Driver et al.|2011| ) with the average error associated with our pho- 
tometric redshifts of the order of Az ^ 0.16. As these photomet- 
ric redshifts apply to both H- ATLAS and non-H- ATLAS sources, 
both populations would experience similar biases associated with 
these errors. Adopting a single redshift at the peak of the photo- 
metric redshift probability distribution (z-PDF) would not accu- 
rately represent our limited knowledge of the redshift of individ- 
ual galaxies. We therefore use the full photometric z-PDF to carry 
out Monte-Carlo (MC) simulations which sample each z-PDF 1000 
times generating 1000 MC cubes for each object. From these sam- 
ples we ensure that we have a good statistical representation of the 
3 -dimensional distribution of galaxies within the survey area. How- 
ever, where available, we use spectroscopic redshifts due to their 
smaller uncertainties. 

3.1 Voronoi Tessellations 

In order to calculate the environmental density of individual galax- 
ies we apply a numerical algorithm called 'Voronoi Tessellations' 
( VT; 'Icke & van de Weygaert 1987 ; van de Weygaert & Icke'1989) 
to the Optical-9hr catalogue. The algorithm works by initially treat- 
ing each object in the field as a single point source (or nucleus). It 
then constructs a convex polygon (or 'Voronoi cell') around each 
of these nuclei enclosing all points that are closer to that nucleus 
than any other. The area of the Voronoi cell is a good representa- 
tion of the local environment of that object, such that the reciprocal 
of this area gives a direct measure of the density. The VT technique 




152 134 136 138 140 142 
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Figure 1. An example of how the VT algorithm works, showing a subset of 
objects placed within the redshift sHce 0.09 < z < 0.10. The VT cells are 
plotted illustrating where the densest regions lie. As the algorithm requires 
more than one object within a specific region to accurately construct a cell 
boundary, the lack of neighbours around the edge of the image results in 
overly large cell areas. As explained in the main text (Section [3T} , these 
outer cells are removed from the rest of the analysis. 



has been used in many areas of astronomy, initially in the study of 
large-scale structure of the universe (e. g.,|Icke & va n de Wey gaert] 
p^87l fvande Weygaert & Icke 1989| [Diehl & Statler 2006) and 
more recently in studies of cluster detection (e.g., Kim et al. 2000; 
Ramella et al.||200T] [van Breukelen et aLl[20b6a, ,van Breukelen] 
et al.|2006bl|van Breukelen & Clewley|2009| [Geach et al.|20lTbl 
Soares-Santos et al.|2 0ir). 

The VT algorithm does not take redshift into account when 
linking each nucleus together as we lack the necessary resolution 
in redshift to apply a three-dimensional VT algorithm. Therefore 
it is necessary to group galaxies into redshift slices so as to avoid 
projection effects. Subsequently each of the 1000 3D MC fields 
are split into redshift slices so that the VT algorithm can be ap- 
plied to each slice individually; thus only objects within each slice 
have the VT algorithm applied to them, maximising the accuracy 
of the density calculation for each object. The width of each slice 
is limited such that, if too wide, projection effects may become a 
problem. In addition, if the slice is too small, an overdense region 
in terms of the rest of the field may go undetected by the VT al- 
gorithm if it is spread across multiple slices. The typical velocity 
dispersion between gravitationally bound group/cluster galaxies is 
within the region of a few hundred kilometres per second ( |Haynes| 
|et al.| 1984] [Martini et al.|2007] ), which equates to a spread in red- 
shift ofAz/(l+z)~10~ atz = l. Therefore we adopt a redshift 
slice of width Az — 0.01, easily incorporating associated galactic 
environments and resulting in 120 slices across each of the 1000 
3D MC realisations across the full redshift range of our data out to 
z 1.2. 

For each object in each MC cube realisation its Voronoi cell 
area {xi) is calculated, the mean of which {xi) gives the overall 
mean area calculated for that object across all of the MC realisa- 
tions. Taking an inverse of this mean area gives a value for the 
mean environmental density {pi) for that object. Figure [T] shows an 
example of one VT slice (containing a smaller subset of the data 
for illustrative purposes) and it is immediately clear that objects to- 
wards the outside of the field have overly large Voronoi cell areas. 
This edge effect is the result of the VT algorithm not finding any 
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objects outside of this boundary and consequently being unable to 
triangulate in these areas. In order to prevent this edge effect alter- 
ing the mean density result, a cut is then made around the outside 
of the field to remove the outermost objects (and their overly large 
cell areas) from the Optical-9hr catalogue. The position of this cut 
was calculated by plotting the right ascension (RA) and declina- 
tion (Dec) values separately against a value that represents a nor- 
malised value for the mean area, the significance {S) (this value is 
introduced to account for a peak in the number density of objects 
as explained in Section [T2| ). The resultant plots showed, unsurpris- 
ingly, a sharp increase in mean area of the cells towards the outside 
of the field and that this edge effect penetrated the field by approxi- 
mately d=0.30 degrees in both RA and Dec. As the Optical-9hr field 
extends well beyond the H-ATLAS-SDP field on all sides, this cut 
is not significant in terms of the number of sources lost and does 
not interfere with the accuracy of our analysis. 

Furthermore, to ensure that the accuracy of the comparison 
between the two samples is maintained, it is necessary to include 
only Optical-9hr objects that reside within the boundary of the H- 
ATLAS-SDP field. This ensures that all objects included in the 
density measure are from across the same region and thus have 
been observed by both SDSS and H-ATLAS observations. Thus 
when comparing far-infrared detected and undetected galaxies we 
are not counting any far-infrared luminous galaxies that would oth- 
erwise be detected in the H-ATLAS SDP catalogue if it were not 
for the boundary limits of the H-ATLAS SDP region. After these 
region cuts are applied, the final Optical-9hr catalogue is reduced 
to 129,518 objects. In section [4?7] this catalogue is divided accord- 
ing to whether or not the galaxies have far-infrared emission in or- 
der that these sub-samples can then be compared. Table [T] shows 
the number of galaxies within these sub-samples in addition to the 
number of galaxies with photometric and spectroscopic redshifts. 



3.2 Density normalization 

With a value for the mean environmental density calculated for 
each object in the Optical-9hr field it is possible to examine the 
3D density distribution across the entire redshift range. Due to the 
flux-density limit of the observations and the much larger volume 
being sampled at higher redshift, a peak in the density distribution 
is found dXz ^ 0.4 corresponding to the peak in the galaxy number 
density. This peak in the detection rate would naturally lead to an 
increase in the mean density being returned by the VT algorithm 
for redshift slices in this range. Therefore two Voronoi cells from 
two different redshift slices cannot be accurately compared in terms 
of their environment. In order to counteract this bias it is necessary 
to normalize the Voronoi cell areas across the entire redshift range 
to produce a normalized environmental density for each object. 

We therefore create a separate random field by applying a ran- 
dom position to each Optical-9hr object within the H-ATLAS SDP 
region within each redshift slice of that MC realisation. We apply 
our VT algorithm to the random field and determine the mean cell 
area xsuce (defined as the sum of the individual cell areas in that 
slice divided by their number) and the standard deviation a of each 
random redshift slice. Our density measure is therefore given by, 

Q Xslice Xi 

— , (ij 

where Xi is the measured VT cell area for each object from the real 
field per MC realisation and Sc is the normalized density value in 
comparison to a random distribution for each object. 




-2 F A 

-2 -1 1 _ 2 

VT Normalised Comparison Density (Sj 

Figure 2. The relationship between the normalized densities returned by the 
VT and NN techniques. The Spearman's Rank Correlation Coefficient (ts) 
between both density distributions returns a value of = 0.961 at a > 5cr 
level, indicating a strong correlation. 

Sc therefore shows how the real field compares to a com- 
pletely random distribution, in terms of the standard deviation a of 
that random distribution, and thus accounts for differences in uni- 
formity of the field between slices. This normalisation also allows 
for the comparison of different objects from across redshift slices 
with different population densities. Taking the mean of these values 
across all MC realisations gives the mean normalized comparison 
density {Sc)- 

3.3 Nearest-neighbour density comparison 

Using VT to probe galaxy environments on individual galaxy scales 
is a relatively new approach to the study of galaxy environmen- 
tal density. Previous work in the analysis of galactic environments 
has instead predominantly relied on estimating the local density 
of galaxies using the projected A^th-nearest neighbour technique 
(Sat), which measures the environmental density in terms of the 
number of galaxies within a circular region defined by the radius 
to the A^th-nearest galaxy (e.g.. Dressier 1980, Lewis et al. 2 002[ 
Mmer et al. 2003, Balogh et al. 2004, Cooper et al. 2005, Silverl 
man et al.,2009, Cucciati et al. 2010, Her nandez-Fernandez et al.| 
2012||Wijesinghe et al.|2012] ). We therefore test our Voronoi Tes- 
sellation density measure (VT) against this A^th-nearest neighbour 
technique (NN) in order to establish whether there are any signifi- 
cant differences between the results obtained from both techniques, 
both in terms of how our overall density correlations are affected as 
well as a comparison of the techniques ability to probe detailed 
structure. 

For our comparison we use A/^ = 5 in line with the majority 
of recent studies that have used the NN technique to examine local 
environmental densities of galaxies (e.g.,|Cucciati et al.|2010[ |Wi-| 
[jesinghe et al.|2012l and |Hernandez-Fernandez et al.|2012| ). We use 
the NN algorithm in exactly the same way as our VT method de- 
scribed in Section [TT] with the NN algorithm applied to each red- 
shift slice within each MC cube, once more normalising the field 
to account for differences in number density and uniformity across 
the redshift range. The only difference between the methods comes 
as a result of the fact that, unlike a VT cell, the S n parameter rep- 
resents larger densities with larger values and thus, to reflect this, 
we reverse the sign of Sc such that positive values once more repre- 
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Right Ascension (deg.) 

Figure 3. Colour co-ordinated plots of the field with all redshift sHces com- 
piled to show the whole density distribution. Red and orange colours rep- 
resent the most overdense regions (positive Sc values) and blue and purple 
represents the most underdense regions (negative Sc values) with the range 
of the normalized density limited to —1 < Sc < 1 for clarity. Top: The 
VT method density output. Bottom: 5th-nearest neighbour method density 
output. The plots confirm that the VT and NN methods reproduce the same 
density structure across the field. 

sent a normalized overdensity. We maintain the boundaries applied 
to our VT density calculation to both prevent such edge effects and 
to maintain comparison accuracy between the two methods. 

We apply both the Voronoi Tessellation and A^th-nearest 
neighbour techniques independently to our Optical-9hr catalogue; 
Figure [2] shows the relationship between the initial outputs returned 
by both techniques. Calculating the Spearman's Rank Correlation 
Coefficient (rs) between both density distributions returns a value 
of rs = 0.961 at a > 5cr level, indicating a strong (although non- 
linear) correlation. In addition, Figure[3]shows the normalized com- 
parison densities (Sc) returned by each technique plotted according 
to RA and Dec positions and coloured according to density. Both 
figures clearly show that each density measure has successfully re- 
produced the same general density structure across the field, with 
the most extreme over- and under-densities located by both meth- 
ods. However, there are some noticeable differences between the 
methods indicating that the intensity of the local density in specific 
regions differs between each method. 

From the direct comparison between the two methods in Fig- 
ure |2] it is clear that the NN method has a greater dynamic range 
in the densest environments where the VT method saturates. Con- 



Table 1. The number of objects within the Optical and FIR sub-samples of 
the initial 129, 515 objects of the Optical-9hr catalogue. These sub-samples 
are also divided according to the number of objects with photometric or 
spectroscopic redshifts in the density analysis. 



Sample 


Total Number 


Number of Photo-z 


Number of Spec-z 


FIR 


2,265 


1,489 


776 


Optical 


127,250 


123,730 


3,520 



versely the VT method appears more suited to distinguishing be- 
tween less dense environments where the NN method saturates. 
Figure[3]shows that using the NN method results in larger regions of 
peak overdensities with less defined regions of intermediate density 
in comparison with the structure distribution from the VT method. 
A full investigation of the pros and cons of different density mea- 
sures has been conducted by |Muldrew et al.|f2012| ) and we refer the 
reader to that paper for more information. But to summarise they 
find that the NN technique is very poorly correlated with the re- 
spective dark-matter-halo mass, although the NN technique is able 
to describe the internal densities of high-mass haloes. 

It is also clear that the initial value selected for N will deter- 
mine the accuracy of the NN technique in various environments. 
Where the value of N remains below the number of associated 
satellites, the measured density will increase with increasing val- 
ues of N. Subsequently in our comparison with the VT method, the 
peak over densities returned from the 5th-nearest neighbour would 
be reduced if, for example, only the 3rd-nearest neighbour were 
used. However, with a larger value of N the NN method will lose 
resolution and become more susceptible to the projected separa- 
tions between distinct overdense regions, influencing the density 
result. 

In contrast, the VT method does not suffer from these issues, 
as essentially the number of neighbours used to define the density 
are not fixed. From the methodology of using the VT algorithm 
(Section [3T| ) it is evident that one does not need to necessarily cat- 
egorise each galaxy into a group or cluster, but can instead simply 
measure the surface density of that galaxy directly from the proper- 
ties of its Voronoi cell. Consequently the VT method is fully adapt- 
able to changes in uniformity of the field and calculates densities 
on individual galaxy scales. Therefore the VT method represents a 
reliable and accurate alternative to the more well established NN 
density measure. 



4 ANALYSIS 

4.1 Far-infrared and control samples 

As described in Section [2?7] we use the likelihood-ratio technique 
of |Smith et al.| ( |201 1\ to associate the far-infrared sources with their 
optical counterparts. This cross-matched sample is hereafter named 
'FIR' (consisting of 2,265 objects) while simultaneously remov- 
ing them from the Optical-9hr catalogue reducing this sample to 
127, 250 objects (hereafter named 'Optical'). These sub-samples 
are shown in Table [T] where they are also divided according to the 
number of galaxies with photometric and spectroscopic redshifts. 

In order to accurately compare how Sc values differ between 
the FIR and Optical catalogues it is necessary to ensure that we are 
comparing like with like, such that the objects selected for compar- 
ison should be considered to be from the same population. By se- 
lecting a matched sample of galaxies based on their colour, SDSS 
r— band model apparent magnitude and redshift distributions we 
ensure that these properties have no influence on any differences 
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Figure 4. g-r vs r-i colour distribution for the 'matched' catalogues Optical 
{red) and Herschel {blue) numbering 2, 706 and 902 sources respectively. 



in environmental density found between the two catalogues. This 
is achieved by gridding the field in four dimensions in order to in- 
corporate diW g — r, r — i, rrir and z parameter space, selecting 
matched objects as only those which share an associated grid space 
in all four planes. The choice of ^ — r and r — i colours are selected 
as we limit our Optical catalogue in the r-band apparent magnitude. 
As explained in Section [22] and as shown in [Taylor et al.| ( [20Tl] ), 
colours provide a reasonable method of matching sources in terms 
of their stellar mass over the redshift range under investigation here, 
although we also investigate this further in Section [43] The grid el- 
ements applied to the total g — r and r — i colour distributions 
incorporate 0.1 and 0.06 magnitudes respectively in colour space. 
This difference reflects the larger range of the total g — r colour dis- 
tribution. Simultaneously, the z and rrir ranges have grid elements 
incorporating 0.02 in redshift and 0.38 in r— band magnitude. 

All Optical sources that share an associated grid space with a 
FIR detection in all four planes are initially grouped as potential 
matches to those FIR objects. Then, within each grid space, a mul- 
tiple of the potential matches (totalling three times the FIR sources 
in that grid space) are selected as matched objects. Selecting Op- 
tical matches equal to a multiple of the FIR sources in each grid 
element allows for a more robust comparison without sacrificing 
any similarities between their distributions. Any additional Optical 
or FIR sources that are not matched are discarded. The FIR sample 
contains considerably fewer objects than the Optical sample (2, 265 
against 127, 250), therefore a large proportion of the Optical sam- 
ple will not have an associated FIR object and thus will be lost 
from the final cross-match. This reduces the sample size to 2, 706 
and 902 for the Optical and FIR samples respectively. 

In addition, the number-density of galaxies reduces with in- 
creasing redshift to such an extent as to affect our sampling. There- 
fore we apply a maximum peak redshift limit onto the samples of 
z < 0.5. This maximum redshift does not influence the Sc values 
of the remaining sample due to the fact that each Sc value already 
incorporates the high redshift galaxies via the full z-PDF sampling 
achieved within the algorithm. Figures [4] & [5] show the colour, red- 
shift and magnitude distributions for these two matched samples. 



4.2 Comparison of the FIR and Optical samples 

We use one- and two-dimensional Kolmogorov-Smirnoff tests (KS- 
tests) to confirm that our matched samples are consistent with hav- 



Figure 5. Redshift vs r— band apparent magnitude (mr) for the 'matched' 
Optical {red) and Herschel {blue) catalogues. Only the redshift range of 
< 2; < 0.5 and rrir range of 15 < rrir < 21.5 was included in the sam- 
pling, outside of these ranges the completeness of the catalogues reduced 
significantly. 



ing been drawn from the same underlying distribution in terms 
of their colour, magnitude and redshift distributions and that the 
various combinations are consistent with each other. These tests 
demonstrate that our null hypothesis, such that the 'Optical' control 
sample is drawn from the same underlying distribution as the FIR 
sample, cannot be rejected at a significant level. The results of these 
tests are presented in Table[2] We then applied the one-dimensional 
KS test to the environmental density measurements for the Optical 
and FIR samples (Sc), which return a probability of just 4.2 x 10~^, 
rejecting the null hypothesis of them being drawn from the same 
underlying distribution at the 3.5cr level. Therefore we find a sig- 
nificant difference in the distribution of the galaxy environmental 
density between far-infrared selected galaxies and a control sample 
with no detectable far-infrared emission. To test this result further 
we applied a two-dimensional KS-test to both Optical and FIR pop- 
ulations comparing the environmental densities in conjunction with 
the optical properties and the redshifts. The two-dimensional KS- 
test comparing the colours, r— band magnitude and redshift distri- 
butions to the environmental density values are presented in Table[2] 
and show that the environmental densities of the FIR and Optical 
samples are significantly different in all cases. 

The extent of this difference is illustrated in Figure [6] where 
we show normalized histograms of the two environmental density 
distributions. It is clear that the Optical population (red), with the 
mean of its distribution at Sc = (12.70 ± 0.95) x 10"^ (denoted 
by the red dashed line), is more overdense (has larger values of Sc) 
than the FIR population (blue), with the mean of its distribution at 
Sc = (3.46 ± 1.49) X 10~^ (denoted by the blue dashed line). Us- 
ing the Mann- Whitney U (MWU) test we also test for differences 
between the median values of the distributions. For the Optical and 
FIR populations the test returns a probability of 5 x 10 indi- 
cating that the two populations have significantly different median 
values at the 4.5(7 level. In Section |431 we show the same result is 
returned when the NN method is used to calculate the environmen- 
tal density, thus showing consistency with our VT result. 
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Figure 6. Normalized histograms showing how the distribution of environmental density (Sc) of the total colour matched Optical (red) and FIR (blue) 
populations compare, with shaded histograms representing the number of spectroscopic redshifts in each sample. Left: Full matched catalogues containing 
2, 706 Optical and 902 FIR objects. The histograms show error bars depicting the normalized error associated with each bin, where Sc > signifies an 
overdensity and Sc < signifies an underdensity in terms of the entire redshift range (0 < z < 0.5). The FIR data are shifted generally to lower Sc values, 
with the mean of the distribution at (3.46 ± 1.49) x 10~^ (blue dashed line), this is contrasted against the mean of the Optical distribution at (12.70 ± 0.95) 
X 10~^ (red dashed line). KS and MWU-tests indicate a significant difference to the 3.5a level. Centre and Right: Normalized histograms that show the full 
matched sample binned in redshift (0 < z < 0.25) and (0.25 < z < 0.50) respectively, showing a continued separation between the distributions increasing 
with redshift from 2.2cr to 3.3cr significance. 



Table 2. Two sample and two-dimensional KS and MWU-test results over 
the full SFR and redshift range (0 < z < 0.5). Where op represents Opti- 
cal (2, 706 objects) and FIR represents FIR (902 objects). The two density 
distributions are significantly different at the 3.5 cr level from KS tests, with 
the medians of the distributions different at the 4.5 cr level from MWU tests. 



Distributions Compared 


KS Prob. 


MWU Prob. 


Zop VS ZpjJi 


0.999 


0.403 


(g - r)op vs (g - r)FiR 


0.974 


0.374 


{r - i)op vs (r - i)FlR 


0.395 


0.239 


^_r{op) vs rnr{FIR) 


0.719 


0.290 


{Sc)op vs {Sc)fir 


< 10-3 


< 10-5 


(g-r,r - i)op ys(g -r,r - i)FlR 


0.224 




{g - r, z)op vs (g - r, z)fir 


0.956 




(r - i,z)op vs (r - i,z)FiR 


0.468 




(rrir, z)op vs (mr, z)fir 


0.755 




{g - r, mr)op vs {g - r, mr)FlR 


0.839 




(r - i, mr)op vs (r - i, mr)FlR 


0.339 




(g - r, Sc)op vs (g - r, Sc)fir 


0.005 




(r - i, Sc)op vs (r - i, Sc)fir 


0.005 




{rrir, Sc)op vs (mr , Sc) FIR 


0.006 




(z, Sc)op vs {z, Sc)fir 


0.010 





Table 3. Full SFR range KS-test and MWU-test results for the comparison 
of both the Optical and FIR populations Sc distributions within individual 
redshift slices. From KS results the density distributions are different be- 
tween the 2.2(7-3.3(7 level. The number of objects from each population are 
given: 

Redshift Slice Optical IR KS Prob. MWU Prob. 

< z < 0.25 680 227 0.028 0.028 

0.25 <z< 0.50 2,026 675 < 10-^ < 10-^ 



4.3 Redshift binning 

Clearly in any flux-density limited sample there is a bias in the 
sense that the higher-redshift sources are more luminous than those 
at lower redshift. Therefore, to further examine the difference be- 
tween the Optical and FIR environmental density distributions we 
split the two populations into two redshift slices ofO < z < 0.25 
and 0.25 < z < 0.5. The < z < 0.25 bin contains 680 objects 
from the Optical population with 227 objects from the FIR popula- 
tion and the 0.25 < z < 0.5 bin contains 2, 026 objects from the 



Optical population and 675 objects from the FIR population. The 
results of the KS -tests comparing the density measurements within 
these bins are shown in Table[3] these show that the null hypothesis 
can be rejected and the two populations can be considered differ- 
ent in terms of their overall density distributions at the 2.2cr level 
for the low-redshift bin and 3.3cr level for the high-redshift bin. 
The MWU test returns probabilities of 0.028 and 3.4 x 10~^ for 
the low- and high-redshift bins respectively, also indicating that the 
two distributions have significantly different median values. These 
binned distributions are shown in Figure [6] Consistent results are 
found when the NN method is used to calculate the environmental 
density, as shown in Section |T6] 

These results show that, as with the full redshift range, ob- 
jects in both redshift bins are significantly different in terms of 
their overall density distributions and their median values. How- 
ever, it is evident that these statistical differences are higher in the 
higher-redshift bin and that this bin contains a larger number of ob- 
jects in both samples. We investigate the impact of this difference 
in number density with increasing redshift by matching the number 
of galaxies in the higher and lower bins and repeating the sample 
comparisons. First, by increasing the redshift boundary between the 
higher and lower bins (from z = 0.25 to 2; = 0.32) to achieve ap- 
proximate matching in sample sizes above and below this redshift. 
Second, by reducing the number of objects within the higher red- 
shift bin to match exactly with the lower bin samples. In both cases 
we find that the same trends are found and our results remain the 
same. 

With fewer objects in the lower redshift bin, the signal to noise 
will be lower at these redshifts, flattening the density distributions 
and affecting the comparison. Furthermore, at higher redshifts, ob- 
jects with lower IR luminosities are excluded by the flux-density 
limit of the H- ATLAS survey. Thus, the density distribution of the 
less luminous far-infrared galaxies may actually be similar to the 
Optical population. In contrast the higher redshift bins contain a 
much higher proportion galaxies with higher levels of star forma- 
tion and therefore exhibit a stronger correlation with density. The 
consequence of this bias is that the statistical differences found be- 
tween the total Optical and FIR distributions (shown in Table [2| 
may be being diluted by galaxies with low SFRs at low redshift. In 
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Table 4. Two sample KS and MWU-test results for each SFR bin, collectively over the full redshift range (0 < z < 0.5) where op represents the Optical 
and FIR represents the FIR populations. (A): SFR of to 15 Mq yr~^ contains 414 cross-matched Optical objects and 138 cross-matched FIR objects. (B): 
SFR of 15 to 30 M0 yr~^ bin contains 1, 104 cross-matched Optical objects and 368 cross-matched FIR objects. (C): Minimum SFR of 30 Mq yr~^ bin 
containing 879 cross-matched Optical objects and 293 cross-matched FIR objects. 



Distributions Compared 


KS Prob. (A) 


MWU Prob. (A) 


KS Prob. (B) 


MWU Prob. (B) 


KS Prob. (C) 


MWU Prob. (C) 


Zop VS ZpjR 


0.999 


0.484 


0.911 


0.359 


0.324 


0.372 


(g - r)op VS (g - r)FiR 


0.834 


0.357 


0.964 


0.420 


0.997 


0.420 


{r - i)op VS (r - i)FlR 


0.054 


0.057 


0.998 


0.453 


0.999 


0.474 


^_r{op) VS rnr{FIR) 


0.994 


0.410 


0.994 


0.368 


0.630 


0.310 


{Sc)op VS {Sc)fir 


0.009 


0.009 


0.006 


< 10-3 


< 10-3 


< 10-5 
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Figure 7. Top: The calculated SFR ( Mq yr~^) against redshift for the total 
FIR catalogue. Bottom: The three SFR bins from the FIR catalogue versus 
redshift. The < SFR < ISM© yr-^ bin (Solid line) containing 781 
objects, the 15 < SFR < 30 yr-^ bin (Dashed line) containing 631 



objects and the SFR 
objects z < 0.5. 



> 3OM0yr-i bin (Dotted line) containing 423 



order to examine the full impact of these objects, in Section|44] we 
apply SFR bins to the FIR catalogue and repeat the above analysis. 

4.4 Star-formation rate vs environmental density 

The SFR of a galaxy can be estimated using the relation given in 
|Kennicutt|(T998 ) as proportional to the total IR luminosity (Lfir) 
over 8-1000//m, assuming a Salpeter IMF between 0.1 M© — 
100 Mq ( |Salpeter|1955t , such that: 



The thermal emission of far-IR galaxies can be represented by a 
modified black body emission spectrum from ,Blain et al.,^2003j : 



,,3+/3 



(If - 1) ^ 



(3) 



SFR{MQyi- 



4.5 X 10" 



I (ergs • 



(2) 



where h is the Planck constant, k is the Boltzmann constant, and 
T represents the temperature. The emissivity index, /3 modifies the 
Planck function by assuming that the dust emissivity varies as a 
power law with frequency, z/^, where /3 can be between 1-2 as de- 
scribed in |Hildebrand| ( |1983| ), depending on the frequency of the 
observations. We fix the dust emissivity index to /3 = 1.5 with dust 
temperature (T) equal to 26 K as found by |Dye et alT] ( |2010| l and 
[Jarvis et aT] ( |2010| ) and integrate the modified blackbody over the 
wavelength range 8 — lOOO/xm to obtain the far-infrared luminos- 
ity. We then use equation|2]to determine the SFR (in Mq yr~^) for 
each galaxy (see Figure [7]r 

As a check of our Lfir values we compare the galaxies in 
our sample to the subset of objects with far-infrared luminosities 
determined from the full energy balance models of da Cunha et aL| 
(2008 ) by Smith et aL| ([2012a). We find that our Lfir values are 
slightly underestimated, and a correction factor of 1.25 is needed 
to produce a 1 : 1 correlation. This suggests that, as we are assuming 
a fixed dust temperature of 26 K, our estimate misses ^ 25 per cent 
of the total dust luminosity, and hence the SFRs are underestimated 
in our calculation. We therefore apply this correction factor to our 
Lfir values to account for this difference in our resultant SFRs. 
We note that this correction makes very little difference to our over- 
all results on the relative environmental densities between different 
bins in SFR. 

It is also worth noting that the SFRsestimated from FIR emis- 
sion may be overestimates of the true SFR. We know that far- 
infrared emission is a tracer of star formation in an idealised case 
where young stars dominate the radiation field and dust opacity is 
high (Kennicutt et al. 2009 ). Multi-temperature dust distributions 
and emission from dust heated by older ISM stars ( da Cunha et al.| 
[2008 1 [Smith et al.|2012a| ) are not be expected to be consistent with 
equation [2] and this also plays a part in our correction factor of 
1.25. This must also be balanced against the fact that any unob- 
scured component of star formation would not be detected in the 
FIR. Thus, although far-infrared emission is highly correlated with 
SFR we note that the absolute values of SFR should be used with 
caution. 

We bin the FIR objects in terms of their SFR in bins of — 15, 
15 — 30 and > SOMoyr"^ in order to compare the impact of 
different SFRs on our initial results. Again matching the control 
sample to the individual binned SFR samples we perform KS- 
tests and MWU-tests on all combinations. The results of these 
tests are shown in Table |4] From KS and MWU tests the density 
distributions for all our SFR bins are shown to be different; for 
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Figure 8. Histograms of the r — K colours for the Optical (red) and FIR 
(blue) samples showing that they have significantly different distributions 
at a > 5cr level with the FIR galaxies lying redward of the Optical galaxies. 



Table 5. Two sample and two-dimensional KS and MWU-test results over 
the full SFR and redshift range (0 < z < 0.5). Where op represents Optical 
(3, 624 objects) and FIR represents FIR (1, 208 objects) matched in terms 
of their band magnitude, r — K colour and redshift distribution. The 
two density distributions are different at the > 5cr level from KS tests, with 
the medians of the distributions different at the > 5cr level from MWU 
tests. 



Distributions Compared 


KS Prob. 


MWU Prob. 


Zop VS ZpiR 


0.996 


0.489 


(r - K)op VS (r - K)fir 


0.626 


0.485 


'^_K{op) VS rnK{FIR) 


0.136 


0.137 


{Sc)op VS (Sc)fir 


< 10-9 


< 10-10 


(r - K, z)op VS (r - K, z)fir 


0.493 




{mK,z)op VS {rriK, z)fir 


0.177 




(r - K, rnK)op vs (r - K,_mK)FlR 


0.313 




{r - K, Sc)op VS (r - K, Sc)fir 


< 10-5 




(mK,Sc)op VS {mK,Sc)FIR 


< 10-'^ 




(z, Sc)op VS (z, Sc)fir 


< 10-5 
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Figure 9. Normalised histograms of the environmental density (Sc) of both 
Optical (red) and FIR (blue) populations, cross-matched inr — K, rriK and 
z parameter space. KS and MWU tests indicate a significant difference to 
the > 5cr level. The FIR data is shifted generally to lower Sc values, with 
the mean of its distribution at 0.09 ± 0.02 (blue dashed line), contrasted 
against the mean of the Optical distribution at 0.23 ± 0.01 (red dashed line) 
Shaded histograms represent the number of spectroscopic redshifts in each 
sample. 

SFR> — 15 M0 yr~^ the difference in environmental density is 
present at the 2.6cr level for both the KS-test and the MWU test, 
for 15 - 30 M© yr~^ the KS and MWU tests shows a difference 
at the 2.7cr and 3.8cr level, and for the SFR> 30 M© yr~^ the dis- 
tributions are significantly different at the 3.3cr and 4.8cr levels, 
respectively. 

This shows that the SFRs derived from the far-infrared emis- 
sion of galaxies are strongly linked to the environmental density 
of the galaxy. Selecting only those galaxies with the highest SFRs 
results in an even more pronounced difference between the normal- 
ized density distributions of the Optical and FIR populations. 

4.5 Dust reddening effects 

In our analysis we have only used the most sensitive optical bands 
to define our optical and far-infrared selected samples, due to the 



wealth of such data over the survey region used. However, we 
would expect that the galaxies which are detected in H- ATLAS to 
be subject to dust reddening effects, as we know that they must have 
significant amounts of dust in them to be detected in the first place. 
This could cause our Optical and FIR galaxies to be mismatched 
in terms of their intrinsic stellar colours and their stellar masses. 
To investigate this we include the i^— band data from the UKIRT 
Infrared Sky Survey (Lawrence et aL|2007|), which is available for 
many (but not all) of our sources. We did not use this initially as 
the number of galaxies with a band identification is less than 
the number which are identified in the ^, r and i bands. In Fig- 
ure [8] we show a histogram of the r — colours of our matched 
sample of Optical (red) and FIR (blue) galaxies where, due to the 
smaller number of i^-band detections, the size of each sample is re- 
duced to 2,139 and 759 objects respectively. This shows that there 
is indeed a significant difference in the r — K distributions be- 
tween the Optical and FIR galaxies at a > 5cr level, suggesting 
that dust reddening may well be biasing our results. However, this 
would only strengthen our results due to the fact that, as we cross- 
match in the r — band, the reddening we see in Figure|8]is caused by 
the FIR population having brighter iC— band magnitudes than the 
Optical sample. Therefore, due to this reddening, we are likely to 
be overestimating the FIR i^— band magnitudes and subsequently 
their masses. Given that we know more massive galaxies generally 
trace denser environments, correcting for this would lead to a larger 
difference between the Optical and FIR samples. 

To test the rigour of our result we repeat our density analysis 
with Optical and FIR samples cross-matched in terms of r—K, ttik 
and z parameter space. As this new cross-matching takes into con- 
sideration only three dimensions, the resultant number of objects 
considered matched in all three of these parameters is larger than 
in our initial four-dimensional cross-matching, with 3, 624 Optical 
and 1, 208 FIR objects. Applying KS and MWU tests to the data 
return probability values consistent with a significant difference be- 
tween the environmental density distributions to the > 5cr level, 
with the FIR population once more favouring underdense regions 
compared to the Optical sample with mean values of 0.09 ± 0.02 
and 0.23 =b 0.01 respectively. Table|5]gives the results of the statis- 
tical comparison and the density distributions are plotted in Figure 

HI 

We do not extend on this analysis here as new data from the 
VISTA VIKING Survey (e.g. [Findlay et aL|2012| ) over the full H- 
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Table 6. Two sample and two-dimensional KS and MWU-test results from 
the application of the 5th-nearest neighbour technique to the SDSS and H- 
ATLAS SDP data from Section[4| Where op represents Optical (2, 706 ob- 
jects) and FIR represents FIR (902 objects). The two density distributions 
are significantly different to the 3.9 cr level from KS tests in agreement with 
our VT technique. 



Distributions Compared 


Jvo rrOD. 


M W U rrOD. 


Zop VS ZpiR 


0.999 


0.403 


(g - r)op vs (g - r)FiR 


0.974 


0.374 


{r - i)op vs (r - i)FlR 


0.395 


0.239 


^r{op) vs rnr{FIR) 


0.719 


0.290 


{Sc)op vs {Sc)fir 


< 10-4 


< 10-5 


(g-r,r - i)op vs(g -r,r - i)FlR 


0.223 




{g - r, z)op vs {g - r, z)fir 


0.956 




(r - i,z)op vs (r - i,z)FiR 


0.468 




(mr, z)op vs (mr, z)fir 


0.755 




{g - r, mr)op vs {g - r, mr)FlR 


0.839 




(r - i, mr)op vs (r - i, mr)FlR 


0.339 




(g - r, Sc)op vs (g - r, Sc)fir 


0.003 




(r - i, Sc)op vs (r - i, Sc)fir 


0.003 




{rrir, Sc)op vs (mr, Sc)fir 


0.003 




(z, Sc)op vs (z, Sc)fir 


0.004 





ATLAS fields will mean that the analysis presented here could be 
carried out with a iT— band selected sample in the near future. 



4.6 Comparison with nearest-neighbour 

Applying the same KS and MWU statistical comparisons from Sec- 
tion |42] between our cross-matched FIR and Optical data sets, we 
find good agreement between the NN and the VT method; the two 
NN defined normalized density distributions (Sc) are significantly 
different, with a KS test probability of 8.7 x 10~^ indicating a 
significant difference at the 3.9a level. As with our VT method 
all other parameter comparisons show no significant difference as 
shown in Table [6] In further agreement with the results established 
using the VT method, the mean values of the two density distribu- 
tions reveal that the cross-matched FIR catalogue contains objects 
with lower environmental densities than the Optical catalogue with 
mean values of 0.31 ± 0.04 and 0.60 ± 0.05 respectively. MWU 
tests reveal a significant difference between the median values at 
the 4.6(7 level. In addition, repeating the redshift binning from Sec- 
tion |43]returns consistent results such that differences are found in 
both bins increasing from 2a in the lowest bin to 3.3cr in the higher 
bin. 



5 COMPARISON WITH SEMI-ANALYTIC MODELS 



In this section we use the semi-analytic models (SAMs) of Hen- 
[riques et ar] ( |2012| ), who construct 24 pencil-beam synthetic light- 
cones for square areas (1.4 deg x 1.4 deg) out to high redshift. 
These light-cones are based on the SAM of |Guo et"aL| p011|) 
which itself is built upon previous SAMs, such as [Croton et at] 
([2006'), to take into account a full range of astrophysical processes; 
reionisation, cooling, disk size, star formation, supemovae feed- 
back, satellites, gas stripping, mergers, bulge formation, black hole 
growth, feedback from active galactic nuclei (AGN), metal enrich- 
ment, dust extinction, stellar mass and the luminosity function. 
Each pencil-beam light-cone is constructed to account for the fact 
that each side of the co-moving Millennium Simulation box is equal 
to 500/i~^Mpc which is smaller than the co-moving distance of 



an object at 2; ^ 2. Therefore periodic replication of the simula- 
tion can lead to multiples of the same object being included due 
to discontinuities in large scale structure at the boundaries between 
replications { [Kitzbichler & White|2007t . 



5.1 SAM analysis 

Here we apply our environmental density measurement to the mock 
catalogues of |Henriques et al.|p012| ) in order to establish whether 
the same relationship between environmental density and SFR is 
found. Our results from Section|4]have shown there is a statistically 
significant difference between the density distributions of galaxies 
with and without obscured star formation (as traced by far-IR emis- 
sion). This is such that galaxies with obscured star formation favour 
underdense regions while galaxies without star formation favour 
overdense regions. 

We use all 24 mock catalogues from|Henriques et al. (2012i, 
with additional data taken from the mock catalogue of Guo et alT] 
pOTT] ). The data are selected to span the same redshift range as 
our observed data (0 < z < 0.5), however these catalogues do 
not initially contain any errors on their redshift values. Therefore 
in order to achieve similar redshift sampling as within our environ- 
mental density measurement in Section[3] it is necessary to apply a 
redshift error to each object based on matching the r— band mag- 
nitude (rrir) and photometric redshift (z) values to the observed 
data catalogue. This therefore establishes a likely value for the red- 
shift error based on these parameters. This is achieved by cross- 
matching each Henriques & Guo (hereafter HG) object to the total 
Optical-9hr catalogue, locating all matches in r— band magnitude 
(rrir) and photometric redshift (z). For each match, the photometric 
error from the Optical-9hr catalogue is applied as the redshift error 
to the HG object. Where a spectroscopic redshift is located, a stan- 
dard error for a spectroscopic redshift (0.001) is applied, as in the 
initial analysis (Section[3]). In addition, as with the observed data in 
Section an apparent r— band magnitude cut was implemented 
removing all galaxies with magnitudes fainter than 21.5. The re- 
sultant number of objects across the HG catalogue totals 260, 303 
objects. Applying our environmental density measure returns nor- 
malized environmental density values in comparison to a random 
field (Sc) for each object, as with our observed data in Section 3.2 



However, due to the smaller field size (—0.7 <RA< 0.7 degrees 
and —0.7 <Dec< 0.7 degrees) the edge effect cuts imposed dur- 
ing the analysis have a greater impact on the number of sources cut 
from our HG catalogue, further reducing the catalogue to 112, 125 
objects. 

For our analysis it is first necessary to determine which galax- 
ies in the simulated catalogue would have far-IR emission, and thus 
be detectable by the H- ATLAS survey. This is achieved by calcu- 
lating the 250/im flux for each object from the given SFR and red- 
shift values given by the SAM. This is the direct reverse of the 
calculation in Section |4.4| where we calculate the SFR, assuming 
a temperature and emissivity index, from the 250/im flux of each 
FIR galaxy. The average 5 a 250/im flux limit of the H- ATLAS 
observations (33.5 mJy), taken from Rigby et al. (2011), provides 
an exact cut-off point for which a galaxy could be considered de- 
tectable in the survey. From here the HG catalogue could be split 
into IR and non-IR detected objects, equivalent to our FIR and Op- 
tical catalogues in our observed data analysis from Section [4T7] Ob- 
jects with a 250 /xm flux density greater than 33.5 mJy are there- 
fore considered detectable by H- ATLAS, hereafter called FIR-HG 
(1, 919 objects), and those with a 250 /xm flux density less than 
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Figure 10. A set of normalized histograms form the SAM output that show how the distribution of environmental density (Sc) of the colour matched Optical- 
HG (red) and FIR-HG (blue) populations change with redshift. Left: The full SFR range of both environmental density (Sc) distributions (2, 184 and 728 
objects respectively), along with error bars depicting the normalized error associated with each bin. The two distributions are significantly different to the 4a 
level as determined by KS tests. Centre: The lower redshift bin (0 < z < 0.22) KS and MWU statistical tests reveal that the distributions are statistically 
different and that the null hypothesis can be rejected to at least the lAcr level. Right: In the higher redshift bin (0.22 < z < 0.47) the distributions are again 
significantly different and the null hypothesis can be rejected to at least the 3.4cr level. 
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Figure 11. Normalized histograms showing the comparison between the 
redshift distributions of both FIR and FIR-HG samples. Solid line: The to- 
tal FIR-HG redshift distribution showing a clear weighting towards lower 
redshifts. Dashed line: The total FIR redshift distribution which peaks at 
higher redshifts. 



Table 7. Full SFR range. Two sample and two-dimensional KS and MWU- 
test results over full redshift range (0 < z < 0.5) where op represents 
Optical-HG (2, 184 objects) and FIR represents FIR-HG (728 objects). The 
density distributions are significantly different to the 4cr level from KS tests, 
with the median values of the distributions different to the 4.1a level from 
MWU tests: 



Distributions Compared 


KS Prob. 


MWU Prob. 


Zop VS Zfjr 


0.658 


0.318 


(g - r)op VS (g - r)FiR 


0.990 


0.452 


{r - i)op VS (r - i)FlR 


0.198 


0.172 


'^_r{op) VS rn^(FlR) 


0.237 


0.125 


{Sc)op VS (Sc)fir 


< 10-4 


< 10-5 


(g -r,r- i)op \s(g-r,r - i)FlR 


0.131 




(g - r, z)op VS (g - r, z)fir 


0.751 




(r - i, z)op VS (r - i, z)fir 


0.278 




(rur, z)op VS (mr, z)fir 


0.081 




(g - r, mr)op vs (g - r, mr)FlR 


0.343 




(r - i, mr)op VS (r - i, mr)FlR 


0.090 




(g - r, Sc)op VS (g - r, Sc)fir 


< 10-3 




(r - i, Sc)op VS (r - i, Sc)fir 


< 10-3 




(rur, Sc)op VS (rur, Sc)fir 


0.002 




(z, Sc)op VS (z, Sc)fir 


0.001 





33.5 mJy are not considered detectable by H- ATLAS and hereafter 
named Optical-HG (110, 206 objects). 

As with our observed data in Section |4.1| a like with like 
cross-matching process is applied. We find 2, 184 objects from the 
Optical-HG population matched with 728 objects from the FIR-HG 
population. Both of these cross-matched samples represent approx- 
imately the same percentage of their parent populations as found 
with our observed cross-matched samples from Section |4.1| This 
was such that the cross-matched Optical and Optical-HG samples 
represent ^ 2 per cent of their parent populations, with the FIR and 
FIR-HG samples representing ^ 40 and ~ 38 per cent respectively. 

5.2 SAM statistical testing 

We perform the same statistical analysis, as with our observed data 
in Section [431 i-^- ^PPlyi^^g one- and two-dimensional KS tests as 
well as MWU-tests to the two cross-matched populations, the re- 
sults of which are given in TablejT] These values show, in agreement 
with our observed data in Section |42] that the null hypothesis can 
be rejected to at least the 4a level for the normalized environmental 



density (Sc) of the Optical-HG and FIR-HG populations. All other 
parameter comparisons cannot be considered significantly different 
at any reliable statistical level (e.g. > 2a). 

The two density distributions are shown in Figure [To|(/^/0. As 
expected, it is the FIR-HG population (blue), with the mean of its 
distribution at = (-4.84±1.54) xlO~'^ (blue dashed line), that 
is biased towards underdense regions, while the Optical-HG popu- 
lation (red), with the mean of its distribution at Sc = (2.59 ±1.72) 
X 10~^ (red dashed line), shows a bias towards overdense regions. 
Errors included on each bin are small and support the difference 
found between the distributions. It is worth noting that if we do 
not include simulated photometric errors and treat the SAM red- 
shift values as precise, the resultant Sc distributions exhibit a much 
larger spread in values. Inclusion of photometric errors evidently 
reduces this spread by essentially 'washing out' the density struc- 
ture we are trying to recover. However, when precise redshifts val- 
ues are used we find the same correlations are found between the 
Optical-HG and FIR-HG populations, despite the increased spread, 
with both populations exhibiting a significant difference. 

To further examine the difference found between the Optical- 
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Table 9. Two sample KS and MWU-test results where op represents Optical-HG and FIR represents FIR-HG. (A): SFR of to 5 Mq yr~^ containing 
Optical-HG 729 objects and 243 FIR-HG objects. The difference between the two density distributions in this bin cannot be distinguished. (B): SFR of 5 to 
10 M0 yr~^ containing 984 Optical-HG objects and 328 FIR-HG objects. The density distributions are different at the 3.3(7 level from KS tests and 3.4(7 
from MWU tests. (C): SFR of > 10 yr~^ containing 669 Optical-HG objects and 223 FIR-HG objects. From KS tests the two density distributions are 
different at the 3(7 level, with the median values of the distributions different at the 3.6(7 level from MWU tests: 



Distributions Compared 


KS Prob. (A) 


MWU Prob. (A) 


KS Prob. (B) 


MWU Prob. (B) 


KS Prob. (C) 


MWU Prob. (C) 


Zop VS ZpiR 


0.175 


0.201 


0.310 


0.277 


0.344 


0.443 


(g - r)op VS (g - r)FiR 


0.947 


0.328 


0.965 


0.423 


0.770 


0.322 


{r - i)op VS (r - i)FlR 


0.125 


0.090 


0.057 


0.080 


0.958 


0.360 


^_r{op) VS rnr{FIR) 


0.065 


0.134 


0.909 


0.301 


0.216 


0.089 


{Sc)op VS {Sc)fir 


0.450 


0.098 


0.001 


0.001 


0.004 


< 10-3 



Table 8. Full SFR range KS-test and MWU-test results for both the Optical- 
HG and FIR-HG populations for Sc distributions within the individual red- 
shift slices shown in Figure [To| The null hypothesis is rejected at the 2.4(7- 
2.9(7 level for both distributions in the lower redshift bin from KS and MWU 
tests. The null hypothesis is rejected at the 3.4(7-3.9(7 level in the higher red- 
shift bin: 



Redshift Slice 


Optical 


IR 


KS Prob. 


MWU Prob. 


0<z< 0.22 


1,273 


418 


0.020 


0.003 


0.22 <z< 0.47 


911 


310 


< 10-3 


< 10-4 



HG and FIR-HG populations we once more split the two popu- 
lations into redshift bins, repeating the analysis from Section [43] 
Figure pT] shows the redshift distributions of both the FIR (dashed 
line) and FIR-HG (solid line) data. It is immediately clear that the 
FIR-HG population is primarily weighted towards lower redshifts 
with a mean value of (1.82 it 0.02) x 10~^. Its number density falls 
off beyond z ^ 0.25 reaching a maximum redshift of 2; = 0.47, 
short of the full redshift range of the FIR sample. In comparison the 
FIR population retains an approximately consistent number density 
across its entire redshift range with a mean value of (3.30 di 0.04) 
X 10~^ . Therefore it is necessary to adjust the redshift binning from 
that applied to the FIR data in Section [43] to account for this differ- 
ence. Reducing the boundary between our higher and lower redshift 
bins from z = 0.25 to ^ = 0.22, we bin the Optical-HG and FIR- 
HG populations, re -plotting histograms to represent the data within 
these redshift bins and reapplying KS -tests and MWU tests to the 
data. The redshift bins are plotted into normalized histograms dis- 
played in Figure [To| and the KS and MWU test results from this 
redshift binning are given in Table [8] These statistical tests confirm 
that the same SFR-density trend with redshift is found as with our 
observational data analysis in Section [43] 

This trend is found despite the differences between the redshift 
distributions of both observed and simulated FIR data sets (Fig- 
ure[TT]). We determine that a galaxy constitutes part of the FIR-HG 
population by selecting objects based on their 250/im flux, which is 
calculated from the individual SFR derived from the SAM of lGuol 
jet al.| { |20lT] ). As noted in Section [44] our calculation of the SFR 
from the far-infrared luminosity is subject to a range of assump- 
tions that may or may not be valid. Furthermore, the parameters 
used within the SAM to calculate SFR may also not accurately in- 
corporate all of the physical processes that govern SFR. Detailed 
analyses of these effects are beyond the scope of this paper, there- 
fore we adopt a conservative approach and consider that the results 
from the observations and the SAMs are in qualitative agreement. 

5.3 Applying SFR limits to SAMs 

Here we repeat the same SFR binning analysis from Section [4.4[ 
Applying this stage of the analysis to the populations derived from 



SAMs allows us to further probe the differences between the sim- 
ulated data and the data obtained observationally. The reduction in 
the number of FIR-HG objects between the low and high redshift 
bins is illustrated in Figure[T2] This shows how the number density 
of the FIR-HG population falls off beyond z ^ 0.25, in compari- 
son with our observed data in Figure[7] It is evident that the FIR-HG 
population has a higher fraction of sources at higher redshifts than 
FIR population. Due to the vast majority of the FIR-HG population 
residing below z ^ 0.25 it is necessary to adjust the SFR binning 
parameters, from those applied to the FIR population, to narrower 
SFR ranges in order to achieve three comparative samples of this 
population. We therefore bin the FIR-HG objects in terms of their 
SFR in bins of - 5, 5 - 10 and > 10 M© yr~\ In agreement 
with our results in Section |4.4| we find that with higher levels of 
star formation, the statistical difference between the two popula- 
tions from KS and MWU tests increases. Only in our lowest SFR 
bin (O-5M0yr-^) is no significant difference found between the 
Optical-HG and FIR-HG samples. These values are presented in 
Table[9] 

Figure [12] shows how at higher SFRs the number density of 
objects reduces to such an extent that it makes further analysis un- 
feasible. Therefore it was not possible to introduce SFR bins at 
higher SFRs than 10 M© yr~^ to the analysis. Despite this, some 
key conclusions can be made with regards to the Optical-HG/ FIR- 
HG comparison based on the three SFR bins applied to the data. 
This analysis has shown that, for galaxies with SFRs higher than 
5 Moyr"^, there is a statistically significant difference between 
both of the Sc distributions and that this statistical difference be- 
comes more pronounced in these higher SFR bins as a result of 
removing lower star-forming objects from the comparison. From 
finding no statistical difference between the density distributions of 
the Optical-HG and FIR-HG populations when SFRs are less than 
5 Mq yr~^ to finding a significant difference to at least the 3a level 
in higher SFR bins. 



5.4 Application of NN to SAMs 

We again test the A^th-nearest neighbour against our Voronoi Tes- 
sellation methods. We apply the NN method to our analysis of 
semi-analytic models. Following the same processes from Sections 
|5.1| and [5.2[ we apply our algorithm, changed to incorporate the 
NN method, to the total HG catalogue, again dividing the output 
according to the H- ATLAS flux limit and cross-matching in ^ — r, 
r — i,z and rrir parameter space. This, once more, provides two cat- 
alogues representative of optical and far-infrared (Optical-HG and 
FIR-HG) that can be accurately compared to analyse differences in 
density. 

We find that our results match those of Section [53] with KS 
and MWU test results finding a significant difference between the 
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Figure 12. Top: SFR ( M© yr-^) versus redshift for the total FIR-HG cat- 
alogue. This shows that above 10 Mq yr~-^ the number of objects reduces 
significantly. Bottom: Plot of the three SFR bins from the FIR-HG catalogue 
vs redshift. The < SFR < 5 M© yr-^ bin (Solid line) containing 846 
objects, the 5 < SFR < 10 M© yr-^ bin (Dashed line) containing 547 
objects and the SFR > lOMoyr-^ bin (Dotted line) containing 526 
objects. 



normalized densities of the Optical-HG and FIR-HG populations 
only (Table [T0|. With a KS-test and MWU-test probabilities of 
^ 10~^ indicating a significant difference at the >4cr level in both 
cases. The mean values of each distribution lie at (4.02 ± 0.24) 
X 10"^ for the Optical-HG and (1.60 =b 0.21) x 10"^ for the FIR- 
HG indicating that the Optical-HG population occupy generally 
more overdense regions in agreement with our study using VT. 



6 DISCUSSION 

The increased statistical separation between the Optical and FIR 
density distributions, found with both increasing redshift and in- 
creasing SFR, provides a clue to the role of environment in the evo- 
lution of galaxies over the redshift range (0 < z < 0.5). We find 
a clear segregation in the galaxy environmental density between 
far-infrared-detected sources and those galaxies that are matched 
in terms of optical colour, magnitude and redshift but devoid of de- 
tectable far-infrared emission. Moreover, we find that this segrega- 
tion becomes more pronounced at brighter far-infrared luminosity 
(or SFR) or at higher redshift. Unfortunately our data precludes us 



Table 10. Two sample and two-dimensional KS and MWU-test results from 
the application of the 5th-nearest neighbour technique to our semi-analytic 
model analysis Optical-HG and FIR-HG populations. Where op represents 
Optical-HG (2, 184 objects) and FIR represents FIR-HG (728 objects). The 
two density distributions are significantly different to the 4.5a level from 
KS tests in agreement with our VT technique. 



Distributions Compared 


JVo rrOD. 


iVlW U rrOD. 


Zop VS Zfjr 


0.658 


0.318 


(g - r)op vs (g - r)FiR 


0.990 


0.452 


{r - i)op vs (r - i)FlR 


0.198 


0.172 


T^jiop) vs rnr{FlR) 


0.237 


0.125 


{Sc)op vs (Sc)fir 


< 10-4 


< 10-^ 


(g -r,r- i)op \s(g-r,r - i)FlR 


0.131 




{g - r, z)op vs {g - r, z)fir 


0.751 




(r - i, z)op vs (r - i, z)fir 


0.278 




(mr, z)op vs (mr, z)fir 


0.081 




{g - r, mr)op vs {g - r, mr)FlR 


0.343 




(r - i, mr)op vs (r - i, mr)FlR 


0.090 




{g - r, Sc)op vs {g - r, Sc)fir 


< 10-3 




{r - z, Sc)op vs (r - i, Sc)fir 


< 10-3 




(mr, Sc)op vs (mr, Sc)fir 


0.001 




(z, Sc)op vs (z, Sc)fir 


< 10-3 





from distinguishing between an evolutionary effect and one associ- 
ated with the level of star formation activity. 

It is important to note that the reliability criterion (R > 0.8), 
that we employ from [Smith et"ar| ( |2011| ) to select optical coun- 
terparts to the FIR data, does not present a bias in our results. This 
potential bias is such that in denser regions, with an increased num- 
ber of potential optical counterparts to a FIR object, the reliability 
parameter for that FIR object, as defined in ,Smith et al.| ( |2012a| l, 
would reduce. In other words, in denser regions it potentially be- 
comes more difficult to reliably associate the FIR detection with a 
unique optical source. Therefore the FIR object may be excluded 
leaving only the FIR detections in relatively low-density environ- 
ments. We tested this bias by making a more inclusive cut to the FIR 
sample based on the minimum likelihood-ratio (LR) rather than the 
reliability criterion (R). Our FIR sample was therefore increased to 
include these previously missing sources. Upon repeating the anal- 
ysis we found that we still obtained similarly significant differences 
between the FIR and Optical samples. 

These results support previous studies that also suggest that 
the presence of star formation in a galaxy is negatively correlated 
with the density of its environment (e.g.. Dressier. 1980, .Postman 



poto 



| & Geller|1984l|Dressler et al.J1997| [Poming uez et al.|2001 
'et al.i 20031|Kauffmann et al.|20041|O'Mill et al.|2008[|Lee et al. 



2010). Our analysis has shown that this correlation holds on in- 
dividual galaxy scales, and thus the processes responsible for this 
correlation must have influence at this level as well as on larger 
scales. In addition, our use of far-infrared observations mean our 
results are not affected by uncertainties associated with extinction 
or Ha to SFR conversions. 

However, the exact mechanism responsible for the observed 
reduction of SFR with increase in density remains uncertain. Re- 
cent studies by [Deng et al.| pOTT] ) and [Wijesinghe et aL] { |2012| ) 
suggest that there is no trend with environment when restricting the 
SFR-density comparison to purely star-forming objects. They con- 
clude, therefore, that the observed SFR-density correlation is due 
to the increasing fraction of passive galaxies across the total galaxy 
sample since z ^ l. |Deng et al.ipoil] ) go further and suggest that 
the SFR-density relation is strongly colour dependent, with blue 
galaxies exhibiting a very weak correlation between environment 
and SFR. In contrast, they find red galaxies to exhibit strong cor- 
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relation between environmental density and SFR attributing this to 
the increasing presence of red late-type morphologies. As our anal- 
ysis has focused on the direct comparison of star formation prop- 
erties with the individual environmental densities of each object, 
we have shown that there is a clear difference between the star- 
forming and passive population in our colour-matched samples in 
agreement with earlier work by Gomez et al.| ( [2003j and Welikak 
|inL| ( [2008l ). 

Furthermore, we have carried out the same analysis on SAMs 
where we obtain a similar result, i.e. the environmental density dis- 
tributions from the total simulated Optical-HG and FIR-HG pop- 
ulations were found to be significantly different at the 4cr level. 
Qualitative agreement is also found when we bin in terms of both 
redshift and SFR. 



led Principal Investigator consortia and with important participa- 
tion from NASA. The H- ATLAS website is: 
http://www.h-atlas.org/ 

GAM A is a joint European- Australasian project based around a 
spectroscopic campaign using the Anglo- Australian Telescope. 
The GAMA input catalogue is based on data taken from the 
Sloan Digital Sky Survey and the UKIRT Infrared Deep Sky Sur- 
vey. Complementary imaging of the GAMA regions is being ob- 
tained by a number of independent survey programs including 
GALEX MIS, VST KIDS, VISTA VIKING, WISE, Herschel- 
ATLAS, GMRT and ASKAP providing UV to radio coverage. 
GAMA is funded by the STFC (UK), the ARC (Australia), the 
AAO, and the participating institutions. The GAMA website is: 
|http ://w w w.g ama- survey, org/ 



7 SUMMARY & CONCLUSIONS 

We have compared the environmental and star formation properties 
of two populations of galaxies out to z ^ 0.5. For this analysis we 
have used optical spectroscopy and photometry from the GAMA 
9hr survey (DRl data) and SDSS, with far-infrared observations 
from the H-ATLAS SDR We use Voronoi Tessellations to analyse 
the environmental densities of these galaxies on individual scales 
normalized to account for differences in the population density and 
uniformity across the redshift range due to the flux limit of the sur- 
vey and the increasing volume sampled with increasing redshift. 

The environmental density of the Optical and far-IR cata- 
logues were then compared by initially matching the catalogues 
in multi-dimensional colour, magnitude and redshift space {g — r, 
r — i, TTir, z) selecting a matched population of the Optical sources 
numbering three times that of the far-IR distribution, in order to 
obtain a robust comparison over < 2; < 0.5. Our key results are: 

(i) Objects with far-IR detected emission, and levels of star for- 
mation > 5 M0 yr~^, reside in less dense environments than galax- 
ies not detected at far-infrared wavelengths. 

(ii) The environmental density difference between the two far- 
IR and non-far-IR luminous galaxies also increases with redshift, 
with a 2.2cr difference in the lower bin (0 < z < 0.25) and a 
3.3(7 difference in the higher bin (0.25 < z < 0.50), with the 
far-infrared detected galaxies again residing in less dense environ- 
ments. In relation to this, we find an increasing separation between 
the density distributions with increasing SFR from 2.6a, 2.1a and 
3.3(7 respectively, although we note that we cannot distinguish red- 
shift effects from luminosity effects in our flux-density limited sam- 
ple. 

(iii) We find substantial differences between redshift distribu- 
tions of both our observed and SAM far-infrared samples. This pro- 
vides interesting indications on how recipes for star formation need 
to be modified within SAMs to improve their ability to model the 
observed universe. 

(iv) We also note that VT are a reliable and accurate method of 
calculating the environmental densities for individual galaxies. In- 
deed, the use of VT for this purpose may surpass the NN technique, 
as their improved resolution is able to measure more detailed den- 
sity structure. 
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